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Решена задача определения температурного поля в поперечном сечении тепловой трубы, проходящем через зону подвода теп- 
ла к внешней поверхности трубы, с учетом неоднородности граничных условий на внешнем контуре. Установлено, что отклоне- 
ния температур по угловой координате не превышают 0, 1 К при существенно неоднородном теплообмене по контуру трубы. 
Максимальные отклонения по толщине испарившегося слоя хладагента также не превышают 1 %. Полученные результаты поз- 
воляют сделать вывод о том, что комплекс теплофизических процессов, протекающих в тепловой трубе, обуславливает практи- 
чески однородное температурное поле в любом сечении корпуса трубы по продольной координате даже при существенном из- 
менении интенсивности теплообмена по внешнему контуру. 


Введение 

Совершенствование современной радиоэлект- 
ронной аппаратуры (РЭА) неразрывно связано с ре- 
шением проблемы обеспечения требуемого тепло- 
вого режима проектируемых устройств. Одним из 
эффективных способов обеспечения теплового ре- 
жима современной РЭА является применение теп- 
ловых труб (ТТ), работающих по принципу замкну- 
того испарительно-конденсационного цикла [1-5]. 

Выбор схемы технической реализации блока 
"источник выделения теплоты - соединительный 
элемент - тепловая труба" представляет сложную за- 
дачу из-за необходимости обеспечения как можно 
более близкого к однородному температурного поля 
в поперечном сечении ТТ. При моделировании теп- 
лопередачи, учитывающей влияние условий подво- 
да энергии к внешней поверхности тепловой трубы 
на температурное поле корпуса ТТ, можно оценить 
эффективность ее работы, а также спрогнозировать 
оптимальные схемы передачи теплоты к ТТ от ис- 
точника тепловыделения. В реальных же условиях, 
из-за неоднородности температурного поля, трудно 
обеспечить, равномерную по окружной координате 
интенсивность подвода тепла в зону испарения ТТ. 
Поэтому теоретический анализ закономерностей 
процесса теплопереноса в рассматриваемой системе 
с учетом изменения условий теплообмена по угло- 
вой координате является актуальным. 

Цель данной работы - численный анализ тем- 
пературного поля в поперечном сечении тепловой 
трубы, проходящем через зону подвода тепла к по- 
верхности ТТ, с учетом неоднородности граничных 
условий на внешнем контуре трубы. 

Физическая модель 

Упрощенная схема блока "тепловая труба - сое- 
динительный элемент" представлена на рис. 1. 

Рассматривается процесс переноса тепла в по- 
перечном сечении ТТ, проходящем через зону ис- 
парения и соединительный элемент (рис. 1), через 
который тепло от РЭА подводится к внешней пове- 
рхности трубы. Такое сечение является наиболее 
информативным, потому что для него характерны 


максимальные тепловые потоки в зону испарения. 
При постановке задачи учитываются процессы пе- 
реноса тепла по радиальной и окружной координа- 
там в ТТ за счет теплопроводности, тепловой эф- 
фект испарения на границе раздела жидкой и паро- 
вой фаз хладагента. При постановке задачи рас- 
сматривались различные варианты граничных ус- 
ловий на внешнем контуре трубы. В данной работе 
представлены результаты исследований, получен- 
ные для самых неблагоприятных по режиму подво- 
да энергии к ТТ граничных условий. Принималось, 
что часть поверхности внешнего контура трубы (I) 
(рис. 1) соединена через элемент (II) с границей, 
через которую энергия от блока РЭА передается на 
тепловую трубу. Оставшаяся часть поверхности 
трубы считалась теплоизолированной. В реальных 
системах обеспечения теплового режима, напри- 
мер космических аппаратов типа спутников связи 
[6], отдельные элементы работают в условиях глу- 
бокого вакуума и отсутствия хорошего теплового 
контакта между собой. Поэтому рассматриваемая 
схема является не только реальной по условиям 
теплообмена на внешнем контуре, но и достаточно 
типичной для перспективных систем. 



Рис. 1 . Упрощенная схема области решения: I) тепловая тру- 
ба; II) соединительный элемент 
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При постановке задачи не моделировались про- 
цессы переноса теплоты в осевом направлении по 
паровому каналу, фитилю и корпусу трубы. Интен- 
сивность перетока энергии в этом направлении не 
зависит от распределения температур по окружной 
координате в рассматриваемом сечении ТТ. Поэто- 
му неучет этого фактора не может повлиять на ос- 
новные закономерности теплопереноса по окруж- 
ной координате в сечении г=сошІ. При этом толь- 
ко необходимо учитывать реальный нестационар- 
ный характер изучаемых процессов. 

В работе [5] задача тепломассопереноса решена 
без учета неоднородности граничных условий по 
окружной координате. В данной работе рассматри- 
вается процесс теплопереноса в блоке "соедини- 
тельный элемент - тепловая труба" в цилиндричес- 
кой системе координат в сечении ^=соп8І. Подра- 
зумевается, что продукты испарения хладагента 
уносятся мгновенно из зоны испарения, а давление 
в паровом канале не зависит от координат г и Ѳ. 

Постановка задачи 
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Система уравнений, описывающих теплопере- 
нос в тепловой трубе в системе координат, с осью 
симметрии, совпадающей с осью симметрии ТТ, 
включает уравнения теплопроводности (1, 2) для 
жидкой и паровой фаз хладагента, а также уравне- 
ние теплопроводности (3) для корпуса ТТ. Теплопе- 
ренос в соединительном элементе описан уравне- 
нием теплопроводности (4) в декартовой системе 
координат. Принималось что корпус трубы и сое- 
динительный элемент изготовлены из одного мате- 
риала с высокой теплопроводностью. 

Геометрия области решения приведена на рис. 2. 



Рис. 2. Г еометрия области решения: 1-2) соединительный 
элемент; 3) тепловая труба 
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С -теплоемкость, ДжДкг-К); Ь -размер зоны под- 
вода теплового потока, м; п - нормаль к внешней 
поверхности ТТ; 0 - теплота фазового перехода, 
Дж/кг; - удельная плотность теплового потока, 
Вт/м 2 ; Я - внешний радиус тепловой трубы, м; 
/•-текущий радиус, м; Г-температура, К; Ж- мас- 
совая скорость испарения, кг/(м 2 с); х, у - коорди- 
наты, м; А - расстояние от источника тепла до 
верхней точки корпуса ТТ на оси симметрии, м; 
Ѳ - угловая координата, град; Я - коэффициент 
теплопроводности, Вт/(м К); р - плотность, кг/м 3 ; 
г - время, с; % - толщина испарившегося слоя жид- 
кости, м. 
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Индексы: 0 - начальный; 1 - паровая фаза; 
2 - зона фитиля; 3 - корпус тепловой трубы; 4 - со- 
единительный элемент. 

Массовая скорость испарения хладагента рас- 
считывалась по формуле: 

у- 

^2кК 0 Т гф /м’ 

где А - коэффициент аккомодации; Р - давление, 
Па; 7? 0 - универсальная газовая постоянная, 
ДжДмоль-К); М -молекулярный вес, кг/моль. 
Индексы: г.ф. - граница фаз; н - насыщенный. 
Задача рассматривалась при следующих основ- 
ных допущениях: 

1. теплофизические характеристики в зоне фитиля 
тепловой трубы рассчитываются как эффективные 
с учётом объёмных долей каждой компоненты; 

2. приход жидкого хладагента в зону испарения 
равен расходу парообразного охладителя в этой 
зоне (процесс массообмена стационарен) и ко- 
ординаты границы испарения не смещаются; 

3. теплофизические характеристики (теплопровод- 
ность, теплоемкость) не зависят от температуры; 

4. вся внешняя поверхность области решения теп- 
лоизолирована, за исключением участка х=0, 
0 <у<Ь (источник подвода тепла); 

5. контакты на границах раздела фаз, на границе 
"корпус тепловой трубы - конденсат", а также 
"труба - соединительный элемент" считаются 
идеальными. 

Принятые допущения не накладывают принци- 
пиальных ограничений на общность постановки 
задачи, и отражают достаточно реальный режим 
работы тепловой трубы. 

Метод решения 

Задача решена методом конечных разностей. 
Дифференциальные уравнения в частных произ- 
водных представляли в виде разностных двумер- 
ных уравнений [7]. Переход на новый временной 
слой реализовывался с помощью двух "дробных ша- 
гов" по схеме расщепления [7]. На первом дробном 
шаге в соединительном элементе рассчитывается 
перенос тепла по координате х, а на втором дроб- 
ном шаге - по координате у с использованием од- 


номерных разностных уравнений. Соответственно 
для ТТ по координатам г и Ѳ. 

Систему одномерных разностных уравнений 
решали с помощью метода прогонки по неявной 
четырехточечной разностной схеме [7]. 

Давление насыщенных паров определяли мето- 
дом Риделя-Планка-Миллера [8]. 

Для соединительного элемента использовалась 
прямоугольная разностная сетка. Для ТТ, представ- 
ляющей собой трехслойный цилиндр, использова- 
лась сферическая разностная сетка. На границах 
между областями проводилось сопряжение прямо- 
угольной и сферической сеток. Особенность реша- 
емой задачи состоит в наличии локально-сосредо- 
точенных источников теплопоглощения высокой 
интенсивности в малой по толщине области, соот- 
ветствующей зоне испарения хладагента. Физичес- 
ки процесс испарения моделируется только в од- 
ной точке разностной сетки по г, соответствующей 
достижению условий испарения. Поэтому при вы- 
боре сеточных параметров особое внимание уделя- 
лось контролю условий сходимости итерационного 
процесса. Шаги по временной и пространствен- 
ным координатам выбирались из условий сходи- 
мости итераций. 

При проведении численных расчетов использо- 
вались теплофизические характеристики для кор- 
пуса ТТ и соединительного элемента из сплава 
алюминия и рабочей жидкости (для жидкой и 
паровой фаз) [8, 9]. 

Результаты и обсуждение 

В таблице приведены распределения темпера- 
тур по угловой координате при различных радиусах 
тепловой трубы (на поверхности ТТ и в области 
раздела сред "жидкость-пар") в момент времени 
г=100 с, #=240 Вт/м 2 . 

Как видно из таблицы, температура снижается 
по мере удаления от оси симметрии с ростом угло- 
вой координаты. Это снижение температуры обус- 
ловлено переносом тепла от более нагретого участ- 
ка ТТ к менее нагретому. Но максимальный пере- 
пад температур между верхней и нижней точками 
поверхности ТТ (рис. 1) составил не более 0,1 К, 
что позволяет сделать вывод о несущественном из- 
менении температуры с угловой координатой. 


Таблица. Распределение температуры по угловой координате, К 
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На рис. 3 приведено распределение толщины ис- 
парившегося слоя хладагента по угловой координате. 



Ѳ, град. 

Рис. 3 . Зависимость толщины испарившегося слоя жидкости 
X от угловой координаты Ѳ 

Максимальная величина испарившегося слоя 
хладагента расположена в области наибольшего 
теплоподвода к границе фитиля ТТ. Это обусловле- 
но более интенсивным испарением жидкого хлада- 
гента (более высокой температурой) в верхней час- 
ти трубы по сравнению с нижней частью ТТ. 

Можно отметить, что отклонения температур от 
максимальных значений составляют доли градуса 
при изменении угловой координаты от 0 до п. По- 
этому в реальных системах обеспечения теплового 
режима на базе тепловых труб этот эффект можно 
не учитывать. Данный вывод подтверждает и зави- 
симость толщины испарившегося слоя хладагента 
от угловой координаты (рис. 3). Отклонения от 
максимального значения % не превышают 1 %, что 
в реальных условиях может не учитываться. 

Заключение 

Полученные результаты позволяет сделать вы- 
вод о том, что специфика физических процессов, 
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протекающих в тепловой трубе, определяет прак- 
тически однородное температурное поле в любом 
сечении по і при изменении угловой координаты 
от 0 до л' и при существенно неоднородных гранич- 
ных условиях на внешнем контуре трубы. Установ- 
ленные в результате численного анализа законо- 
мерности обусловлены комплексом совместно 
протекающих физических процессов, основными 
из которых являются поглощение энергии при ис- 
парении хладагента, перетекание тепла по корпусу 
трубы и по фитилю как в радиальном, так и в ок- 
ружном направлениях. При этом неоднородность 
граничных условий на внешнем контуре трубы не 
оказывает практически значимого влияния на тем- 
пературное поле системы в целом и на интенсив- 
ность процессов испарения. 

Подвод тепловой энергии к поверхности внеш- 
него контура трубы может быть локализован на дос- 
таточно малом участке границы; важно, чтобы при 
этом обеспечивался отвод энергии от источника 
тепловыделения с необходимой интенсивностью. 

Установленные в результате теоретического 
анализа закономерности создают объективные 
предпосылки и для вывода о возможности матема- 
тического моделирования процессов тепломассо- 
переноса в тепловых трубах на базе осесимметрич- 
ных моделей [10], учитывающих процессы перено- 
са энергии, массы и количества движения по ради- 
альной и осевой координатам. 

Полученный вывод применим для описания 
процессов, протекающих как в низкотемператур- 
ных, так и в высокотемпературных тепловых трубах. 
Переход на режимы высоких температур не вносит 
принципиальных изменений в физические модели 
процессов и, соответственно, не может привести к 
существенным изменениям в режимах теплопере- 
носа только за счет, например, кратного увеличения 
уровня рабочих температур тепловой трубы. 
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